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Abstract. We present here an analysis of the light curves 
Q ■ of 5.3 million stars in the Small Magellanic Cloud observed 
^3 . by EROS (Experience de Recherche d'Objets Sombres). 
^ ■ One star exhibits a variation that is best interpreted as 
. . , due to gravitational microlensing by an unseen object. 
^ ■ This candidate was also reported by the MACHO collab- 
^ , oration. Once corrected for blending, the Einstein radius 
' crossing time is 123 days, corresponding to lensing by a 
?3 i Halo object of 2.61 8 ,' 2 M©. The maximum magnification 
is a factor of 2.6. The light curve also displays a peri- 
odic modulation with a 2.5% amplitude and a period of 
5.1 days. Parallax analysis of the candidate indicates that 
a Halo lens would need to have a mass of at least 0.3 M©, 
although a lens in the SMC could have a mass as low as 
0.07 Mq. We estimate the optical depth for microlensing 
towards the SMC due to this event to be ~ 3.3 x 1CU 7 , with 
an uncertainty dominated by Poisson statistics. We show 
that this optical depth corresponds to about half that ex- 
pected for a spherical isothermal Galactic Halo comprised 
solely of such objects, and that it is consistent with SMC 
self-lensing if the SMC is elongated along the line-of-sight 
by at least 5 kpc. 



Key words: : Galaxy: halo, kinematics and dynamics, 
stellar content - Cosmology: dark matter, gravitational 
lensing 



1. Introduction 



Send offprint requests to: Nathalie.Delabrouille@cea.fr 

* Based on observations made at the European Southern Ob- 
servatory, La Silla, Chile. 



Ten years after Paczynski's proposal ( Paczynski 198€ ) to 
use gravitational microlensing as a tool for discovering 
dark stars, and four years after the identification of the 
first candidate events in the direction of the Large Mag- 
ellanic Cloud (LMC) (Alco ck et al. 1993, Aub ourg et al. 
1993) and Galactic Bulge ( Udalski et al. 1992 ), searches 
for microlensing events have started to yield quantitative 
information that contributes to a better understanding of 
Galactic structure (Stanek et al. 1996). Probably the most 
intriguing result is that the measured optical depth for mi- 
crolensing towards the LMC implies a total Galactic Halo 
mass in compact objects that is within a factor of two of 
that required to explain the rotation curves of spiral galax- 
ies (Alcock et al. 1997c , see also Ansari et al. 1996a). The 
time scales associated with these events indicate surpris- 
ingly high mass lenses and the difficulty of accounting for 
the events with known stellar populations has stimulated 
interest in star formation and evolution processes. Strong 
limits have been set on the maximum contribution of low 
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mass objects to the Halo of the Milky Way (Renault et al. 
1997, see also Alcock et al. 1996). 

Given the importance of these results, it is impera- 
tive to verify them by using other lines of sight, the most 
promising ones being the Small Magellanic Cloud (SMC) 
and M31. Here, we present a first analysis of microlensing 
data in the direction of the SMC by using 5.3 million light 
curves collected by EROS2 during the first year of the sur- 
vey. More details can be found in (Palanque-Delabrouille 
1997). 

2. Experimental setup 

Our results have been obtained with a completely re- 
designed setup. The EROS program now uses exclusively 
the dedicated 1 meter MARLY telescope, specially refur- 
bished and fully automated for the EROS2 survey (Bauer 
et al. 1997), now in operation at the European Southern 
Observatory at La Silla, Chile. The telescope optics al- 
lows simultaneous imaging in "blue" (A in 420 — 720 nm, 
peak at A ~ 560 nm) and "red" (A in 620—920 nm, peak at 
A ~ 760 nm) wide pass-bands of a one-square-degree field. 
This is achieved by a beam-splitting dichroic cube with a 
CCD camera mounted behind each channel. Each camera 
contains a mosaic of 8 Loral 2048 x 2048 thick CCD's. 
The total field is 0.7 deg (right ascension) x 1.4 deg (dec- 
lination). The pixel size is 0.6 arcsec, and typical global 
image quality (atmospheric seeing + instrument) is 2 arc- 
sec FWHM. 

The read-out of the entire mosaic is done in parallel, 
controlled by Digital Signal Processors, and takes 50 sec- 
onds. The data are first transferred to two VME crates 
(one per color), which manage the real-time part of the 
acquisition system, and then to two Alpha workstations 
where a quality assessment is run (monitoring CCD de- 
fects, sky background, seeing, number of stars. . . ) and 
flat-field reduction is done. The raw and reduced data are 
finally saved on DLT tapes. 

Data taking with the new apparatus began in July, 
1996. Microlensing targets include fields near the galactic 
center, in the disk of the Galaxy, and in the LMC and 
SMC. The data discussed here concern 10 fields covering 
the densest 10 deg 2 of the SMC, as illustrated in figure 
|. The fields were observed from July, 1996 to February, 
1997 and again starting in July, 1997. During the 1996-97 
season, from 60 to 120 usable images were taken of each 
field, giving a sampling time of one point every 2-4 days on 
average. Exposure times varied from 5 min in the central 
fields to 15 min in the outermost fields. 

The DLT tapes produced in Chile are shipped to 
the CCPN (IN2P3 computing center, CNRS) in Lyons, 
France, where data processing occurs. For each of the 
fields, a template image is first constructed by adding to- 
gether 10 exposures of good quality, each re-sampled by 
a factor of 0.7. A reference star catalog is then built us- 
ing the CORRFIND star finding algorithm (the stellar de- 
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Fig. 1. Position of the SMC fields (smOOl to smOlO) on the 
de Vaucouleurs' isophote map (de Vaucouleurs 1957). The 
star in field 5 indicates the position of the microlensing 
candidate. The isophote levels (inner to outer ones) are 
21.4, 21.7, 22.0, 22.5, 22.8, 23.2, 23.9 mag/arcsec 2 in red 
light (A = 625nm). 

tection is done on a correlation image, obtained as de- 
scribed in section 3, page||). For each subsequent image, 
after geometrical alignment to the template catalog, each 
star identified on the reference catalog is fitted together 
with neighboring stars, using a PSF determined on bright 
isolated stars and imposing the position from the refer- 
ence catalog. A relative photometric alignment is then 
performed, assuming most stars do not vary. Photomet- 
ric errors are computed for each measurement, assuming 
again that most stars are stable, and parameterized as a 
function of star brightness and image sequence number. 
Photometric accuracy is in the 8 — 20% range at magni- 
tude V ~ 20 (depending on image quality), and of the 
order of 2% for bright stars (V ~ 17) .[] The number of 
reconstructed stars varies from 8 x 10 5 deg -2 in the dens- 
est region (where errors are dominated by crowding) to 
4 x 10 5 deg -2 in the outer regions (where errors are domi- 
nated by signal- to- noise) . The ph otometry is described in 
more details in ( [Ansari R. 1996b[ ). 



3. Data analysis 

The 5.3 million light curves are subjected to a series of se- 
lection criteria and rejection cuts (globally called "cuts") 
to isolate microlensing candidates (Palanque-Delabrouille 
1997). The first three (1-3) make use of the expected gen- 
eral characteristics of microlensing candidates: single vari- 



1 This analysis does not require absolute calibration and the 
magnitudes mentioned here are only approximate. Absolute 
calibration is in progress. 
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ations on otherwise constant light curves, which coincide 
in time for data taken in both colors. The next two cuts 
(4 and 5) are designed to reject a known background of 
variable stars, while the last two (6 and 7) improve the 
signal-to-noise of the set of selected candidates. The cri- 
teria were sufficiently loose not to reject events affected 
by blending or by the finite size of the source, or events 
involving multiple lenses or sources. We define a positive 
(negative) fluctuation as a series of data points that (i) 
starts by one point deviating by at least la from the base 
flux, (ii) stops with at least three consecutive points below 
base flux + la (above base flux — la) and (iii) contains at 
least 4 points above base flux + la (below base flux —la). 
The significance LP of a given variation is defined as the 
negative of the logarithm of the product, over the data 
points it contains, of the probability that each point devi- 
ates from the base flux by more than the observed fluctu- 
ation (xi is the deviation of the point taken at time ij , in 
er's, N is the number of points within the fluctuation): 

«.--|>(iE*(^)) (1) 

We order the fluctuations along a light curve by decreas- 
ing significance. The cuts of the analysis are described 
hereafter: 

— 1: The main fluctuation detected in the red and blue 
light curves should be both positive and occur simulta- 
neously: if I is the time interval during which the data 
are more than la away from the base flux, we require 

(Ired (1 Ibluc)/(Ircd U Ibluc) > 20%. 

— 2: To reject flat light curves with only statistical fluc- 
tuations, we require that on a given light curve 
LP(2 nd most significant fluct.) / LP(main fluct.) < 
0.35 

in both colors 

— 3: We require that iP(main fluct.) > 30 in both colors. 

— 4: To exclude short period variable stars which exhibit 
scattered light curves, we require that the RMS of the 
distribution of the deviation, in cr's, of each flux mea- 
surement from the linear interpolation between its two 
neighboring data points be smaller than 2.5. 

— 5: We remove two under-populated regions of the 
color-magnitude diagram (see figure ||) that contain 
a large fraction of variable stars ((3 Cephei, RV Tauri 
variables, semi-regular giant variables and Mira Ceti 
stars), defined by (flux F given for the EROS2 filters 

— R for red and B for blue — normalized to an expo- 
sure time of 300 s): 

\og(F R /F B ) < -0.20 and log(F fl ) > 4.5 
\og(F R /F B ) > +0.07 and log(F R ) > 2.7 

2 This does not reject multiple lenses with caustic crossings 
since all the points magnified would be contained in the same 
main fluctuation. 



Fig. 2. Cut on color-magnitude diagram, here shown for 
field smOOl. The dotted lines delimit the rejected areas. 
The dots correspond to all the light curves in the field, the 
star markers are the remaining objects for this field, after 
cuts f through 4. 



— 6: We remove events with low signal-to-noise by re- 
quiring a significant improvement of a microlensing fit 
(ml) over a constant flux fit (est), i.e. that 

[X 2 (cst) - x 2 (ml)]/[x 2 (ml)/d.o.f.] > 150 

where d.o.f. is the number of degrees of freedom. 

— 7: We require that the maximum magnification in the 
microlensing fit be greater than 1.40. 

The tuning of each cut and the estimate of the effi- 
ciency of the analysis is done with Monte Carlo simulated 
light curves. To ensure similar photometric dispersion on 
simulated events and on the data, the events are added to 
real light curves. The microlensing parameters are drawn 
uniformly in the following intervals: time of maximum 
magnification to £ [tfirst — 150, ti ast + 150] days, impact pa- 
rameter normalized to the Einstein radius uq € [0, 2] and 
time-scale (Einstein radius crossing time) At £ [0, 150] 
days. We correct for blending statistically, using a study 
of the typical flux distribution of the source stars which 
contribute to the flux of a reconstructed star, depending 
on its position in the color-magnitude diagram. 

Table | summarizes the impact of the cuts. The first re- 
quirement removes 98% of the data light curves which are 
just flat light curves. It also removes a large fraction of the 
simulated events of too low amplitude (large impact pa- 
rameter) or short-duration events peaking well outside the 
observational period [ifi rs t, iiast]- The other cuts remove a 
large fraction of remaining data light curves (background) 
while leaving, in general, at least 75% of the simulated 
light curves. 
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Cut description 


Number of 

stars 
remaining 


Fraction of remaining 
stars removed by cut 


Data 


Simulation 


Stars analyzed 


5,277,858 






JL . kJllllLllLdllCl IV 


125,071 


98% 


80% 


2: Uniqueness 


36,032 


71% 


13% 


3: Significance 


4,022 


89% 


13% 


4: Stability 


1,214 


70% 


16% 


5: HR diagram 


463 


62% 


6% 


6: Microlensing fit 


48 


89% 


20% 


7: Magnification 


10 


76% 


16% 



Table 1. Impact of each cut on data and simulated events. 
Each fraction for cut n refers to the stars remaining after 
cut (n — 1). 



The efficiency of the analysis (cuts 1 through 7) for 
detecting real microlensing events is determined from the 
set of simulated microlensing events, taking into account 
the effect of blending. The efficiencies (in %) normalized 
to an impact parameter uq < 1 and an observing period 
T bs of one year are summarized in table for various 
Einstein radius crossing times At (in days). 



At 


7 


22 


37 


52 


67 


82 


97 


112 


127 


e(Af) 


8 


16 


20 


22 


24 


27 


28 


29 


30 



At 


150 


300 


500 


1000 


1500 


2000 


e(Af) 


29 


28 


27 


24 


19 


17 



Table 2. Efficiency (in %) of the analysis (cuts 1-7) for 
various time-scales At (in days), normalized to uq < 1 
and T b s = 1 yr. We monitor A ob s = 5.3 x 10 6 stars. The 
efficiency values for At > 150 days are obtained by simu- 
lating events of duration At over a period [tfirst — At, t\ ast + 
At] days. 



4. Study of the candidates 

Of the 5.3 million light curves, ten events passed all cuts 
and were inspected individually. Scanning of Monte Carlo 
events indicates that a negligible number of remaining 
candidates would be rejected by visual inspection. Three 
of the ten candidates exhibit new variations on their light 
curve when adding the first data from the second year 
survey, and are probably recurrent variable stars. Another 
event has its light curve affected by the appearance of a 
neighboring object which is below our detection threshold 
in the template image, but becomes very bright for a 
period of about 60 days. The object could be a nova or 
even a microlensing on an unresolved star. This analysis, 
however, uses only stars identified on the template image, 
so the light curve must be rejected. Four more events 
have light curves incompatible with microlensing and are 
probably due to other p hysical processes (o ne of them 
is a nova in the SMC ( Alcock et al. 1997a )). Another 
event exhibits a very chromatic variation (Ar = 1.7 and 



Ab = 1.2). If this were due to blending (i.e. when both 
the magnified star and an undetected star contribute to 
the total reconstructed flux), the star undergoing the 
magnification would be of similar brightness as clump 
giant stars but redder than clump giants by about 1.2 
magnitude. It could not belong to the SMC and a mi- 
crolensing interpretation of this event is thus unrealistic. 
Thus, 9 out of the 10 candidates are rejected by this 
inspection. 

The remaining light curve is shown in figure H It fits 
well the standard microlensing hypothesis with an Ein- 
stein radius crossing time of 104 days, a maximum mag- 
nification of 2.1 (impact parameter uq — 0.53) occurring 
on January 11, 1997 and x 2 / d -°- f - = 268/161 = 1.7. 



EROS - MICROLENSING CANDIDATE SMC #1 
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Fig. 3. Light curve of microlensing candidate SMC #1, 
with a standard microlensing fit (combined for the red 
and blue light curves) superimposed, with no blending as- 
sumed. Time is in days since Jan. 0, 1990 (Julian date 
2,447,891.5). 

The source star is located at a = 01 h 00'5.64" and 
S = -72°15'2.41" (J2000), and is labeled U0150_00676152 
in the USNO star catalog. The magnitude of the source 
star given in the catalog is R ~ 18. This microlensing 
candidate was also reported by the MACHO collabora- 
tion and exhibited no variation during 3 years preceding 
the upward excursion detected here (Alcock et al. 1997b). 

Because of the high stellar density of the fields mon- 
itored in microlensing surveys, the flux of each recon- 
structed star generally results from the superposition of 
the fluxes of many source stars. We thus introduce two 
additional parameters in the fit : the contribution Cbi in 
both colors of the base flux of the magnified star (/ s tar) 
to the total base flux recovered (/ s tar + /blend): 



Cbl 



/si 



^reconstructed 



1 



/star H~ /blend 



A T 



cal 



(2) 
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The blending coefficient Cbi is unity when there is no 
blending and Cbi — > in the limit where the magnified 
star does not contribute at all to the total recovered base- 
line flux. Allowing for blending does not significantly im- 
prove the fit, but changes the best estimate values of the 
fit parameters, as shown in table [| The errors on the 
blending coefficients are of 20%. The magnified source star 



Fit 


uo 


to 


At 


cw R 


Cbj B 


xVd.o.f. 


Red 


0.53 


2567 


100 






164/94 


Blue 


0.54 


2567 


103 






99/64 


Combined 


0.53 


2567 


104 






268/161 


Blended 


0.41 


2567 


123 


0.71 


0.71 


266/159 



Fig. 4. Template image and correlation image around the 
candidate. North is to the left, East is up. 



the light is split with the ratio 70% to 30% between the 
two blended components. This independent method thus 
gives a result consistent with that of the microlensing fit. 

As can be seen in figure ^, the measured luminosi- 
ties exhibit an abnormally high scatter in the two colors. 
Correspondingly, the x 2 /d.o.f. (266/159 in Table 2) has 
a low probability, of order 10 -5 . As a significant correla- 
tion is observed between the residuals of the fit in both 
colors (99.5% CL), a search for periodicity was performed 
on these residuals. The most likely period was found at 
P = 5.124 days, with a false detection probability of 
2. 10 -6 (other periods, aliases of P, are less probable, 
though not excluded) . Figure || shows the residuals light 
curve, folded to P = 5.124 days. About 1.5% of the main 
sequence stars of similar brightness have a light curve 
which exhibits a periodic modulation of at least a ~2% 
amplitude. 

We then repeated the microlensing fits of table 2, in- 
cluding a sinusoidal modulation with three additional free 
parameters: period, phase and amplitude of the modu- 
lation (identical for both colors). Results of the fit are 
given in table ||. The fits give almost identical results for 
Mo, to, At and blending factors as those in table 2 (respec- 
tively 0.42, 2568., 123. and 0.74). Of course, the modu- 
lation can affect either the amplified component, or the 
blended companion. The results of the fits favor the first 
possibility, though only at the 2.5cr level. We expect that 
the data reported by the MACHO collaboration (Alcock 
et al., 1997b) have enough statistical power to discrimi- 
nate between these two possibilities. We remark that the 
X 2 of the fits including a modulation term are satisfactory, 
indicating accurate modeling of errors. 

Figure |^ illustrates the position of the candidate re- 
constructed star in the color-magnitude diagram of the 
surrounding region (star marker), as well as that of the 



Table 3. Results of microlensing fits to the SMC candi- 
date, to is the time of maximum magnification and At the 
Einstein radius crossing time, both given in days. 

would then have log(fn) = 3.90 and log(/fl//s) = —0.38 
while the blend companion would have log(/^) = 3.53 
and log(/fl//s) = —0.39. This amount of blending is in 
agreement with the estimate given by the MACHO col- 
laboration. 

The blending hypothesis is strengthened by the ob- 
servation of a slight shift in the position of the centroid 
of the reconstructed star as the magnification occurs: 
(AS ~ 0.072 arcsec, Aa ~ -0.090 arcsec). This shift can 
be caused by the displacement of the barycenter of the po- 
sitions (weighted by the flux) of the two components of the 
blend. The direction of this displacement (slope \Aa/AS\ 
of 1.25) is compatible with the direction of the apparent 
elongation of the source star (AS ~ —0.95 arcsec, Aa ~ 
1.26 arcsec i.e. a slope of 1.3) in the image of the correla- 
tion coefficients between a Gaussian Point Spread Func- 
tion (PSF) and the template image (the correlation image 
is shown in figure [|, right): 

_ covariance (PS F, Image) 

coeff. = v 7 == y ' = (3) 

■J variance (PSF) variance (Image) 

The correlation is calculated over a radius of 1.6<7, where 
the a of the PSF is related to the seeing by seeing = 
2v2 In 2 a. The reconstructed source star might therefore 
consist of two components, located approximately 1.5 arc- 
sec apart. Both the template image and the correlation 
image of an 8 arcsec x 8 arcsec area around the candidate 
are shown in figure |. The pixel size on the template image 
is 0.42 arcsec, that on the correlation image is 0.21 arcsec. 
Note the clear improvement in stellar separation (and in 
stellar detection, as evidenced with the leftmost star) in 
the correlation image as compared to the template image, 
which allows us to infer the existence of a blend companion 
and estimate its position. Requiring on the template im- 
age (figure [|, left) the existence of two source stars located 
along the observed position angle (instead of a single star 
recovered as with the standard star finding algorithm), we 
can estimate the flux ratio of the two stars. In the best fit, 
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Fig. 5. Light curve of the residuals in red, folded to P = 
5.124 days. 



Star modulated 


^4mod. 

(in %) 


(in days) 


X 2 /d.o.f. 


magnified star 
blend companion 


2.9 ± 0.5 
11. ± 7. 


5.128 ±0.004 
5.128 ±0.004 


157/156 
163/156 



Table 4. Result of microlensing fit + sinusoidal modula- 
tion on either component of the blend. A mo( j. is the am- 
plitude of the modulation normalized to the unamplified 
stellar flux, and -P m od. the period. 



two components of the blend (circles) . They all lie on the 
main sequence. 

5. Estimate of optical depth and lens mass 

The optical depth is the instantaneous probability that a 
given source star be magnified by more than a factor of 
1.34. It can be estimated as 



1 



N ohB T ohs 2^ ts e(At) 



E 



At 



(4) 



where e(At) is the detection efficiency given in table pi 
T bs = 1 year and iV b s = 5.3 10 6 stars. With the char- 
acteristics of the single event described above, this yields 
(fit with blending): 



3.3 x 10" 



(5) 



i.e. about 50% of the optical depth predicted by a "stan- 
dard" isothermal and isotropic spherical halo fully com- 
prised of compact objects (cf section 5). It is consistent 
with that measured toward the Large Magellanic Cloud 
(Alcock et al. 1997c, see also Ansari et al. 1996a). 

Assuming a standard halo model with a mass frac- 
tion / composed of dark compact objects having a single 
mass M, a likelihood analysis allows us to estimate the 
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Fig. 6. Color-magnitude diagram of the field surround- 
ing the microlensing candidate. The exposure time on this 
field was 480 s and is here renormalized to 300 s for com- 
parison with figure ^. The recovered star is plotted with a 
star marker, the components of the blend with circles. 

most probable mass of the deflector generating the ob- 
served event. The likelihood is the product of the Pois- 
son probability of detecting N evt events when expect- 
ing / N M , by the probability of observing the time-scales 
(Aii, ■■ j Af cvt ). We calculate likelihood contours in the 
(log(M),/) plane using a Bayesian method with a uni- 
form prior probability density in / and in log(M) (i.e. 
equal probability per decade of mass). They are shown in 
figure 0. We integrate over / to obtain the 1-D likelihood 
for the mass of the deflector. This yields the most probable 
mass of the Halo deflector, given with la error bars: 

M = 2.6+«; 2 3 M (6) 

More statistics is obviously required to constrain the mass 
of halo deflectors. 



6. Expected number of Halo events — model 
dependence 

We studied a wide range of disk-halo models. Table || sum- 
marizes their characteristics. The rotation velocity at the 
Sun predicted by these models is always within the obser- 
vational range: VTot(R©) ^ 200km/s (Merrifield 1992), or 
VTot(Ro) = 220 ± 15km/s (Binney and Tremaine 1987). 
The optical depth is independent of the mass of the deflec- 
tors; the event rate is given assuming that all deflectors 
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Fig. 7. Likelihood contours with one microlensing candi- 
date for a standard halo model. The cross marks the peak 
of the 2-D distribution (solution with blending). 

in the halo have a single mass equal to 1 M Q . To obtain 
the predicted value of the event rate for other masses, one 
only needs to scale by t]m, the integral of the mass de- 
pendence of the event rate times /m{M), the normalized 
mass distribution: 

Vm = [ f M (M)-=t== dM (7) 

This simplifies to [Mj M©) -1 / 2 for a Dirac distribution 
peaked at M. 
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V T ot(R Q ) (km/s) 
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205 
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189 
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134 


163 


Wot (50 kpc) 


199 
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180 


148 
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r(10~') 


6.8 


5.7 


2.1 


3.9 


4.2 


3.8 


r(10- 7 yr" 1 ) 


22.8 


17.8 


5.4 


14.2 


12.6 


9.1 



Table 5. Description of the Galaxy models. We give the 
mass M of the halos out to the SMC (in units of 10 11 M©), 
the rotation velocities (in km/s), the optical depth r and 
the event rate T for 1M Q deflectors (with a 100% effi- 
ciency). So is the central column density, /?© the local 
halo density at the Sun, /? is proportional to the asymp- 
totic slope of the rotation curve and q is the flattening 
ratio of the halo. 



Model 1 is the "standard" halo model: an isotropic and 
isothermal spherical halo, with a mass distribution given 
in spherical coordinates by (Caldwell and Coulson 1981): 



where R c = 5 kpc is the Halo "core radius" and i?© = 
8.5 kpc is the distance from the Sun to the Galactic Cen- 
ter. Model 2a is the equivalent power-law model (Evans, 
1993). Model 2b has a maximal disk (E = 100 M /pc 2 ) 
and a very light halo, model 5 intermediate disk and halo. 
Model 3 has a flattened E6 halo (axis ratio q — 0.71) 
and model 4 a decreasing rotation curve (/3 = 0.2 where 
/3 is proportional to the logarithm of the asymptotic 
slope). Models 2-5 are all power-law halo models, with 
self-consistent mass and velocity distributions. For each 
of these, we derive the expected number of events versus 
the mass M of the objects in the Halo (assuming a Dirac 
mass distribution) and compare with the observations (cf 
figure |h. 




p(r) = Po 



R Q 



m 



(8) 



og(M/M„) 



Fig. 8. Expected number of events for the models de- 
scribed in the text, including the experimental efficiency. 



Whatever the halo model considered, the sole event 
we observed, if caused by a deflector in the halo of our 
Galaxy, corresponds to at least 40% of the total optical 
depth expected for / = 1. This is a very small number 
statistics, however. 

7. Parallax analysis 

The very long time-scale of the observed event suggests 
that it could show measurable distortions in its light curve 
due to the motion of the Earth around the Sun, (the par- 
allax effect: Gould 1992), provided that the Einstein ra- 
dius projected onto the plane of the Earth is not much 
larger than the Earth orbital radius. The first detection 
of parallax in a gravitational microlensing event was ob- 
served by Alcock et al. (1995). The natural parameter to 
measure the strength of parallax is the semi major axis 
of the Earth orbit, Ro, in units of the projected Einstein 
radius: Su — i?o(l — x)/Re where x — Dd/D s , with Da 
the distance from the observer to the deflector and D s the 
distance to the source. 

No evidence for distortion due to parallax is detectable 
on the light curve, implying either a very massive deflector 
with a very large Einstein radius, or a deflector near the 
source. Because our results for the standard and blended 
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fits (see table |^) agree well with those obtained by the 
MACHO collaboration with a 3 year baseline (Alcock et 
al. 1997b), we will fix the level of the baseline flux to that 
obtained previously with the blended fit, to perform par- 
allax fits. We fit simultaneously the red and blue light 
curves allowing for parallax and for the periodic modula- 
tion described in section 3. Assuming a blending coefficient 
Cbi = 0.74, our data allows us to exclude, at the 95% CL, 
that Su > 0.054. This yields a lower bound on both the 
projected transverse velocity of the deflector: 



PARALLAX CONSTRAINT ON LENS MASS AND DISTANCE 



5 = Ro/{At Su) = v t /(l -x)> 270 km/s, 
and on the projected Einstein radius: 
R E = Atv = R E I {I - x) > 18 AU. 
We can thus write: 



(Re) 2 



AGM 



x D s 



M 



M P> X 1 



1-x 



> 0.7 



> (18 AU) 2 



(9) 
(10) 

(11) 
(12) 



The high projected transverse velocity definitely excludes 
the possibility that the deflector is in the disk of the 
Milky Way, where the typical velocity dispersion is ~ 
30 — 40km/s (Binney and Tremaine 1987). Moreover, a 
disk lens (i.e. x < 1/100) generating this event would have 
a mass M > 70 M Q ! For a deflector in the halo, x < 2/3 
at the 95% confidence level (for a standard halo) which re- 
quires the mass of the deflector M to be at least 0.3 M , 
while for a deflector in the SMC, if 1 — x ~ 1/10, the mass 
of the lens would be M ~ 0.07 M©. This is illustrated in 
figure ^. 

It is possible for some parallax distortions to be largely 
cancelled out by blending effects. However, blending dis- 
tortions of light curves are always symmetric about the 
point of highest magnification, while this is not the case, 
in general, of parallax distortions. It is only true when 
the velocity of the deflector is parallel or anti-parallel to 
the velocity of the Earth around the Sun at the moment 
of highest magnification. Figure ^ illustrates the amount 
of blending required to compensate the effect of an in- 
creasing parallax while remaining compatible with the ob- 
served light curve. All the points plotted yield a x 2 /d.o.f. 
for the fit within la of the minimum value (157/158). 
Note the two minima regions in the planes shown in the 
figure, one around an angle of 180 degrees between the 
projected velocities of the Earth and of the deflector (full 
markers), while the other (empty markers) corresponds 
to a null angle. The shaded area delimits blending coeffi- 
cients greater than unity, which is not physical. As shown 
in figure [l0| Su could be larger than 0.054, but only in 
the unlikely case of alignment of velocities. We built a 
likelihood based on the fit with parallax, blending and a 
modulation on the magnified star, taking into account the 




0.7 0.8 0.9 1 

D(lens) / D(source) 



Fig. 9. Relation mass-distance of the deflector, from the 
parallax analysis. Only the region above the curve is al- 
lowed. The top gray area is the la most probable mass of 
a Halo deflector (see section 4). 



PARALLAX AND BLENDING 
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Fig. 10. Degeneracy valleys between parallax and blend- 
ing effects. Top plot: blending coefficient compensating the 
effect of parallax. Bottom plot: angle (in degrees) between 
the directions of the line of sight and deflector velocities 
compensating the effect of parallax. Full and empty mark- 
ers correspond to the same x 2 /d.o.f., but the first set tends 
to an angle of 180 degrees, while the second tends to an 
angle of 0. The four sets of markers correspond to the four 
possible relative motions of the Earth and the deflector. 



probability that the alignment of the Earth and deflec- 
tor velocities were parallel or anti-parallel. This yields the 
95% CL upper limit Su < 0.06 (which requires a blending 
coefficient Cbi = 0.45). The projected velocity is then con- 
strained to be at least 190 km/s and equation |lj becomes 
x t% > 0.5. 



M 

M Q 
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8. Discussion — SMC self lensing 

If the deflector belongs to the Halo of our Galaxy, it is 
expected to have a mass greater than a third of a solar 
mass (see sections 5 and 7) ; and yet to be dim enough to 
avoid direct detection, it could only be a white dwarf, a 
neutron star or a black hole. It is also possible, however, 
that both the lensing object and the source star belong to 
the SMC, in which case the deflector would have a much 
smaller mass. 

Let us estimate the optical depth for SMC self-lensing. 
Various authors have suggested that the SMC is quite 
elongated along the line-of-sight, with a depth varying 
from a few kpc (the tidal radius of the SMC is of the order 
of 4 kpc) to as much as 20 kpc, depending on the region 
under study (Hatzidimitriou and Hawkins 1989, Caldwell 
and Coulson 1986, Mathewson et al. 1986). We will ap- 
proximate the SMC density profile by a prolate ellipsoid: 



2h 



£ e -\ z \/ h e - r / r d 



where z is along the line-of-sight and r is transverse to the 
line-of-sight. The depth h will be a free parameter, allowed 
to vary between 2.5 and 7.5 kpc. Fitting for instance the 
Mathewson et al. Cepheid data (Mathewson et al. 1986) 
in the bar of the SMC with the above density distribution 
gives h ~ 5.8 ±1.2 kpc (assuming Poissonian statistics on 
the Cepheid counts). The other parameters are estimated 
from the surface-brightness map of de Vaucouleurs (see 
figure [j]) using the identity: 



R = 26.1 mag/arcsec 2 



1 L Q /pc 2 



(14) 



derived from M v = 4.83 - 2.5 log(Ly/L 0! y), Mr — My = 
—0.35 for the Sun. Plotting the isophote levels R as a 
function of the mean distance of each isophote to the op- 
tical center of the SMC (see figure [ll]) , we can derive the 
central surface brightness: Rq = 20.7 mag/arcsec 2 . The 
slope of the fit yields the value of the radial scale length: 
rd = 0.54 kpc. 

Assuming a mass-to-light ratio of ~ 3Mq/L q , this 
gives a central surface density So — 400 M©pc -2 and a 
total SMC mass of ~ 1 x 10 9 M Q , compatible with that 
estimated from the mass of the LMC considering that the 
SMC is only ~ 20% as bright. We denote as Zd and z s 
the positions of the deflector and the source, both in the 
SMC, with origin taken at the center of the SMC. 

Assuming the same spatial distribution for the source 
stars and the lenses, the optical depth for SMC self-lensing 
is given by 



Jr(z s ) 



M 1 



dr dO dz„ 



I 



hp- r dr dO dz s 



(15) 



Such a large scale-length in the depth of the SMC has been 
criticized, however, by Martin et al. (1990). 




1.25 1.5 
mean radius (kpc) 



Fig. 11. Isophote level (from the de Vaucouleurs' map) 
vs. mean distance to optical center of the SMC. 



where t(z s ) is the usual optical depth due to source stars 
all at a distance z s and deflectors at Zd contained in the 
elementary volume r dr d9 dzd'- 



(13) t(z s ) 



p(z d ) , AGM (D s 

dZd 7T 

M c? 



+ Z d )(z s - Z d ) 



{D s + z s ) 



(16) 



For h — 2.5, 5.0 or 7.5 kpc, this yields optical depths r = 
1.0 10 -7 , 1.7 10~ 7 or 1.8 10 -7 respectively. Considering 
the very limited statistics we have, this optical depth is 
consistent with the observations. 

Let us also consider the expected typical time-scales of 
SMC-SMC microlensing events. The velocity dispersion in 
the SMC is <cr> - 30 km/s (Hatzidimitriou et al. 1997, 
Suntzeff et al. 1986), so the estimated mass M of the de- 
flector causing the observed event (At = 123 days) could 
be greatly reduced compared to that of a halo deflector. 
On average, we have: 



M 

— — x x(l 

M© 



0.0088 



(17) 



Thus, if the deflector is 5 kpc (resp. 2.5 kpc) from the 
source, we have M ~ 0.1 M© (resp. 0.2 M©). 

As more data are accumulated, we expect SMC-SMC 
events to be highly concentrated in high density regions of 
the Cloud (see figure 1), unlike Halo events which should 
be distributed like the SMC stars over the sky. In ad- 
dition, they should not have measurable parallax distor- 
tions. These criteria will help distinguish between the two 
possibilities. 

9. Conclusion 

We have presented here the result of a one-year survey to- 
ward the Small Magellanic Cloud with EROS2. One star 
has a light curve that is best interpreted as due to mi- 
crolensing with an Einstein radius crossing time of 123 
days when allowing for blending, with 70% of the total 
baseline flux contributed by the star being lensed. The 
light curve exhibits a 2.5% modulation with a period 
P ~ 5.12 days. The optical depth estimated from this 
event is ~ 3.3 x 10~ 7 , to be compared with r = 6.8 x 10 -7 
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for a spherical isothermal halo containing only such ob- 
jects. The most probable mass of the deflector (if it is in 



the Halo) would be M — 2.6I2 3 M©. ^ we interpret this 
event as due to a deflector in the SMC, the expected op- 
tical depth is r ~ 1.7 x 10 -7 for a depth scale-length of 
the SMC of 5 kpc, and the mass of the deflector would be 
reduced to about 0.1 M©. Given the very small statistics, 
the SMC interpretation seems possible. Furthermore, we 
observe no distortion on the light curve due to the varying 
velocity of the Earth on its orbit around the Sun (paral- 
lax), although some would have been expected for such a 
long duration event, unless the deflector were either very 
heavy or near the source. This further supports the SMC 
lens interpretation. 

Further observations will help discriminate Halo from 
Cloud deflectors. In particular, because the velocity dis- 
persions in the LMC and the SMC differ by almost a factor 
of 2, the observation of a significant trend for longer time- 
scale events toward the SMC than toward the LMC would 
be a clear signature of events dominated by SMC or LMC 
self-lensing. 
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Note added in proof After submission of this article, the 
OGLE-2 collaboration confirmed, with their data taken after 
June 1997, the modulation we detected. They also separated 
the two components of the blend, which allowed them to es- 
tablish that t he modulation affe cts the magnified component 
(Udalski et al. |astro-ph/9710365| ). The MACHO collaboration 
provided us their data on the candidate, which allowed us to 
confirm the period, P — 5.124±0.001 days, and the amplitude, 
(2.4±0.4)%, of the modulation; we thank them for this private 
communication. 
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